Strange behavior of persistent currents in small Hubbard rings 
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Abstract 

We show exactly that small Hubbard rings exhibit unusual kink-like structures giving anomalous oscilla- 
tions in persistent current. Singular behavior of persistent current disappears in some cases. In half-filled 
systems mobility gradually drops to zero with interaction, while it converges to some finite value in non- 
half-filled cases. 
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Since the pioneering work of Biittiker, Imry, and 
Landauer , the phenomenon of persistent current in 
mesoscopic systems has been discussed theoretically 
in various papers El H 03 HE El El M EDI EH E21 EH 
Although this phenomenon is thought to be qual- 
itatively understood within the framework of one- 
electron theory [Q El El E 03 El j it has proved difficult 
to explain recent experimental results [TH IT51 ITTfl H7] . 
The unexpected large amplitude of the measured cur- 
rents, the observation of currents with both hc/e and 
hc/2e flux periodicity, and, the sign of the currents 
are the major subjects of disagreement. There are 
many theoretical attempts to explain the experimen- 
tal results, and it is generally believed that disor- 
der El El EHl E] and electron-electron correla- 
tion [3 El El EDI EH E1E31 are the key factors which 
can resolve the controversy between theory and ex- 
periment. However, no consensus on the role of inter- 
action and disorder has yet been reached. The expla- 
nation of the experimental results become even much 
more illusive since recently Kravtsov et al. EH have 
shown that additional currents may be generated in 
the rings by other mechanisms which are not experi- 
mentally distinguishable from persistent current. 

To understand the precise role of correlation, in 
this letter we concentrate on the exact calculation 
of persistent current and Drude weight in ordered ID 
Hubbard rings threaded by magnetic flux. Our study 
reveals that certain aspects of the many-body effects 
on persistent current have not been discussed in the 
literature, and we will see that Hubbard correlation 
leads to many remarkable effects. We limited our- 
selves to small Hubbard rings and these studies might 
be also helpful to understand the physical prop- 
erties of TTF-TCMQ conductors, various aromatic 
molecules and systems of connected quantum dots jHJ- 
With the new technological advancements, it is now 
possible to fabricate such small rings, and, in a recent 
experiment Keyser et al. EH reported the evidence 
of anomalous Aharonov-Bohm oscillations from the 
transport measurements on small quantum rings with 
less than ten electrons. In these small rings with few 
number of electrons, the electron-electron interaction 
becomes important as the Coulomb potential is not 
screened much, and we will show that the electron- 
electron interaction leads to similar anomalous oscil- 



lations in the persistent current as a function of the 
magnetic flux. 

We use Hubbard model to represent the system, 
which for a iV-site ring enclosing a magnetic flux <fi 
(in units of elementary flux quantum (f> = hc/e) can 
be expressed as 

H = t H [e Je 4, CT c *+i, CT + &.c] + U n n n U (!) 



where = 2n(f>/N, and the other symbols have their 
usual meanings. We take t = — 1 and use the units 
c = 1, e = 1. 

At zero temperature, the persistent current in the 
ring threaded by flux <p is determined by [3] 



(2) 



where E ((j)) is the ground state energy. We eval- 
uate this quantity exactly to understand unambigu- 
ously the role of interaction on persistent current, 
and this is achieved by exact diagonalization of the 
many-body Hamiltonian Eq. ftp. 

We also consider the response of the system to a 
uniform time-dependent electric field in terms of the 
Drude weight [201 |2"T] D, a closely related param- 
eter that characterizes the conducting nature of the 
system as originally noted by Kohn [221 • T ne Drude 
weight can be calculated from the relation E0] 



D = 



N (d 2 E{ct)) 
4^ 



(3) 



where 4>min is the location of the minimum of Eo(<f>). 
Now we present below the results of our calculations 
for the various systems. 

We first consider two electron systems, and to have 
a deeper insight to the problem, we begin our dis- 
cussion with the simplest possible system which can 
be treated analytically up to certain level. This is 
the case of a three-site ring consisting of two oppo- 
site spin electrons. The Hamiltonian is a (9 x 9) 
matrix which can be block diagonalized to (6 x 6) 
and (3 x 3) matrices by proper choice of the ba- 
sis states. This can be done by constructing ba- 
sis states for each sub-space with a definite value 
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Tabic 1: Eigenvalues (A) and eigenstates of two opposite spin electrons for N=3. 
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of the total spin S. The basis set for the six- 
dimensional sub-space is chosen as A= {| T4,0, >, 

|o,u,o >, |o,o,tl>, (I U,o > -| |,t,o >)/V2, 
(|o,T,|> -|o,|,t>)/72, (I l,o,t> -I T,o,|>)/V2} 

(each with S = 0), while that for the other sub- 
space is taken as B = {(| |,l,0 > +| |,T,0 >)/72, 

(|o,T,|> +|o,|,t>)/72, (| T,o,|> +| |,o,T>)/72} 

(each with S = 1). For U = and U = 2, in Ta- 
ble we list the energy eigenvalues, eigenstates and 
the degeneracy of the levels of this system in the 
absence of any magnetic field. In this table we set 
a = -72/(1 + 73) and 6 = 72/(-l + 73). Both for 
[7 = and U ^ cases, the two-fold degeneracy of 
the energy levels corresponds to the C% v symmetry 
of the system, and, this degeneracy gets lifted by the 
magnetic flux. It is apparent from this table that the 
eigenvalues, eigenstates and also the degeneracy of 
the levels are not affected by correlation in the three- 
dimensional sub-space (the reason is that basis set B 
does not involve any doubly occupied site), but this 
is not true in the other sub-space. 

The insertion of magnetic flux <f> does not alter the 
above structure of the Hamiltonian, though it be- 
comes now field dependent. In Fig. ^a) we have 
plotted some of the low-lying energy levels -E(</>)'s as 
functions of <f> (say, E — 4> curves) of this three-site 
ring with [7 = 0, and, the corresponding I(4>) versus 
4> curve (say, I — cf> curve) is shown in Fig. db) . The 
E — (f) and 7 — curves for U — 2 case are respec- 
tively given in Fig. ^c) and Fig. ^d). In Fig. ^a) 




Figure 1 : The E — <j) and I — cf) curves for following 
two (|, {) electron systems : i) N = 3, U = in (a) 
& (b); ii) N = 3, U = 2 in (c) & (d); and iii) N = 25, 
U = 3 in (c) & (f). 

and Fig. OJc), the dotted curves represent the U- 
independent energy levels as mentioned above. Quite 
interestingly Fig. \%£c) shows that even in the pres- 
ence of interaction, one of the U -independent energy 
levels becomes the ground state energy level in cer- 
tain intervals of <j> (e.g., around <j) — ±0.5). In other 
regions of <j> we have the expected results, i.e., the 
ground state energy increases due to correlation and 
the Eq — cf> curves become much shallow. As a result, 
the usual sawtooth shape of the I — <fi characteris- 
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tics of U = case changes drastically as shown in 
Fig- Old), and the role of correlation is quite com- 
plex rather than a simple suppression of the persis- 
tent current as predicted earlier. A sudden change in 
the direction and magnitude of persistent current oc- 
curs, solely due to correlation, around 4> = ±0.5, and, 
it forms kink-like structures in the persistent current 
as illustrated in Fig. ^d). The kinks become wider 
as we increase U, and, the most surprising result is 
that the persistent current remain invariant inside the 
kinks irrespective of correlation. It is also seen that 
correlation does not affect the <j>$ flux periodicity of 
the persistent current. 

With the above background, we next investigate 
the role of correlation on persistent current in larger 
mesoscopic rings. In Fig.^e) and Fig.[TJf) we display 
respectively the E—(f> and I—<f> curves for a ring con- 
taining two opposite spin electrons with TV = 25 and 
C/ = 2. These figures respectively resembles Fig.^c) 
and Fig.^d), and it implies that correlation plays the 
same role in both the cases. Thus we can precisely 
say that away from half- filling, a ring consisting of 
two opposite spin electrons always exhibits kinks in 
the persistent current for any non-zero value of U. 

Next we consider rings with two up and one down 
spin electrons as representative examples of three 
electron systems. In Fig. EJa) and Fig. HJb), we re- 
spectively plot the E — <p and I — 4> curves for the 
half-filled situation (N = 3 and n = 3, where n 
denotes the number of electrons) with U — 3. We 
find that correlation just diminishes the magnitude 
of the current compared to that of the noninteract- 
ing case. In Fig. |2Jc) and Fig. |2fd) we respectively 
plot the E — 4> and I — <fi curves of a non-half-fillcd 
situation with N = 10 and [7 = 3, while Fig. Efe) 
and Fig. E£f) are these curves for the U = 15 case. 
The behavior of the persistent current, as a function 
of flux, are quite different at low and high values of 
U. It is evident from Fig. EJc) that for low U, all 
the [/-independent energy levels (dotted curves) al- 
ways lie above the ground state energy of the sys- 
tem, and Fig. |2Jd) shows that, apart from a reduc- 
tion of the persistent current, the I — </> curve looks 
very similar to that of a ring without any interac- 
tion. On the other hands, we see from Fig.Ete) that 
for high U, one of the [/-independent energy levels 
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Figure 2: The E — <fi and I — cj> curves for following 
three (T,t,i) electron systems : i) N = 3, U = 3 in 
(a) & (b); ii) AT = 10, U = 3 in (c) & (d); and iii) 
N = 10, U = 15 in (e) & (f). 



becomes the ground state energy in certain intervals 
of 4> (e.g., around <f> — 0), and, this produces kinks 
in the I — <j> characteristics as depicted in Fig. EI 1 )- 
Certainly there exists a critical value U c of correla- 
tion above which kinks appear in the I — cj> curve. 
We have verified that these features of the energy 
spectra and persistent current are the characteristics 
of any non-half-filled ring with two up and one down 
spin interacting electrons. We also note that the half- 
filled cases exhibit </>o/2 periodic currents, while the 
non-half-filled cases have 0o periodicity. 

Now we consider rings with two up and two down 
spin electrons as illustrative examples of four electron 
systems. Let us first study the half-filled case (N = 4 
and n = 4). In the absence of any interaction, this 
system exhibits sharp discontinuity in persistent cur- 
rent at certain values of 4>. However, the influence of 
correlation is quite dramatic and it makes I{4>) a con- 
tinuous function of </> as depicted in Fig.^b). This is 
a quite fascinating result since correlation drastically 
changes the analytic behavior of 1(4)), and we will 
see that this result also holds true for other half-filled 
systems with even number of electrons. Away from 
half-filling, we study an interesting typical case of a 
six-site ring with two up and two down spin electrons. 
It may be considered as a doubly ionized benzene- like 
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Figure 3: The E — (f> and I — <f> curves for following Figure 4: The E — <j) and I — <fr curves for following 
four (T,T,I,I) electron systems : i) N = 4, U = 2 five (T,ti1\ 1,1) electron systems : i) N = 5, U = 2 
in (a) & (b); ii) N = 6, U = in (c) & (d); and hi) in (a) & (b); ii) N = 6, U = 2 in (c) & (d); and hi) 



iV = 6, U = 6 in (e) & (f). 



AT = 6, U = 120 in (e) & (f). 



ring, a system with the promise of experimental veri- 
fication of our predictions. In Fig. Etc) and Fig.|3jd) 
we respectively plot the E — <p and I — <f> curves with 
U = 0, while Fig. I^e) and Fig. OJf ) are these plots 
with {7 = 6. Fig. Off) shows that kinks appear in the 
I—<j) curve (e.g., around = 0) for any non-zero value 
of U, but we interestingly note that the kinks are now 
due to [/-dependent eigenstates. Both the half-filled 
and non-half-filled rings exhibit 0o periodic currents. 

As specific cases of five electron systems, we study 
rings with three up and two down spin electrons. The 
E — (j> and I — cf> curves for the half-filled ring (N = 5 
and n = 5) are plotted in Fig. 0Ja) and Fig. Efb) 
respectively. Interaction just diminishes the magni- 
tude of the current and we find <po/2 periodicity. As 
a non-half-fillcd five electron case, we take a singly 
ionized benzene-like ring (N = 6, n = 5) and calcu- 
late persistent current for low and high values of U. 
In Fig. QJc) and Fig. 0{d) we respectively draw the 
E — (j> and I — 4> curves taking U = 2, while Fig. Etc) 
and Fig. Eff) are these diagrams with U = 120. It is 
clear from Fig. Eld) and Fig. Elf) that like the three 
electron cases, kinks appear in the persistent current 
only after a critical value of U. This reveals that the 
characteristic features of persistent current in meso- 
scopic Hubbard rings with odd number of electrons 



are almost identical. 

Lastly, we consider six electron systems, and inves- 
tigate the nature of persistent current in rings with 
three up and three down spin electrons. At half-filling 
(a benzene-like case with N = 6 and n = 6), we 
see that /(</>) becomes a continuous function of <f> as 
shown in Fig. |5jb), exactly similar to the half-filled 
four electron case. Fig.OJc) and Fig.EJd) are respec- 
tively the plot of E — 4> and I — <f> curves for a typical 
non-half-filled case with N = 7 and n = 6, and, we 




Figure 5: The E — <fi and I — tf> curves for following 
six (|, t, t, J., I, I) electron systems : i) N = 6, U = 2 
in (a) & (b) and ii) N = 7, U = 2 in (c) & (d). 
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find striking similarity in the behavior of persistent 
current with other non-half-fillcd systems containing 
even number of electrons. Hence it becomes apparent 
that mesoscopic Hubbard rings with even number of 
electrons exhibit similar characteristic features in the 
persistent current. 

In the present context it is also interesting to study 
the conducting properties of these mesoscopic rings 
in the presence of Hubbard correlation. For this pur- 
pose, we calculate Drude weight D of these systems 
using Eq. J3J|. A metallic phase (ballistic regime) will 
be characterized by a finite value of D, while it will 
be zero in an insulating phase (localized regime) |22| . 
In Fig. Eta) and Fig.|^b), we plot D as a function of 




ducting irrespective of the value of U. 

In conclusion, we studied exactly the persistent 
current and Drude weight in small mesoscopic Hub- 
bard rings at zero temperature, and, identified that 
the Hubbard correlation leads to many interesting 
new results. The main results are : the appear- 
ance of kinks in the persistent current, observation 
of both tfio and </>o/2 periodic currents, no singular 
behavior of persistent current in the half-filled sys- 
tems with even number of electrons, evidence of U- 
independent eigenstates, existence of both metallic 
and insulating phases, etc. We also observe discon- 
tinuities in persistent current at non-integer values 
of 0o due to electron-electron correlation, which cru- 
cially depends on the filling of the ring and also on the 
parity of the number of electrons. This gives rise to 
anomalous Aharonov-Bohm oscillations in the persis- 
tent current with much reduced period where period- 
icity is not perfect, and Keyser et al. ^H] experimen- 
tally observed similar anomalous Aharonov-Bohm os- 
cillations in the conductance of a few-electron quan- 
tum ring. Finally, it is worth noting that the above 
mentioned effects of interaction will be pronounced in 
strongly correlated electron systems, like, mesoscopic 
rings of transition metal elements. 

We thank R. K. Moitra and B. Bhattacharya for 
many illuminating discussions. 



Figure 6: Drude weight vs. interaction for (a) half- 
filled and (b) non-half-filled systems. The solid, 
dashed and dotted lines respectively corresponds to 
3, 4 and 5 electron cases. 

U respectively for half-filled and non-half-filled rings 
with n = 3,4 and 5. In the non-half-filled cases, the 
number of sites corresponding to a given value of n is 
taken as N = n + 1, and, we note that other choices 
of N do not affect the basic features of the D versus 
U curves. It is evident from Fig.[|Ja) that for low U, 
the half-filled systems are in the metallic phase, and, 
they become insulating only when U is very large. In 
the insulating phase ground state do not support any 
empty site, and the situation is somewhat analogous 
to Mott localization in ID infinite lattices. Fig.^b) 
shows that the non-half-filled rings are always con- 
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